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Abstract 

We present a QCD study on B, D — > ir semileptonic transitions at zero momentum transfer 
and an estimate of magnitudes of the associated CKM matrix elements. Light cone sum rules 
(LCSRs) with chiral correlator are applied to calculate the form factors f§~* n (0) and f?~~* n (0). 
We show that there is no twist-3 and-5 component involved in the light-cone expansions such that 
the resulting sum rules have a good convergence and offer an understanding of these form factors 
at twist-5 level. A detailed 0{a s ) computation is carried out in leading twist-2 approximation 



and the MS masses are employed for the underlying heavy quarks. With the updated inputs 
and experimental data, we have /f ^(0) = 0.281^2 and \V uh \ = (3.4±g;§ ± 0.1 ± 0.1) x 10~ 3 ; 
/f "^(0) = 0.62 ± 0.03 and \V cd \ = 0.244 ± 0.005 ± 0.003 ± 0.008. As a by-product, a numerical 
estimate for the decay constant fn is yielded as fo = 190t;[iMeV. 
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I. INTRODUCTION 



An intensive study on the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements remains 
a cornerstone of high energy physics programme, in testing the standard model (SM) and 
exploring new physics. The unitarity of the CKM matrix must be put to a test by phe- 
nomenological research on the so called unitarity triangle. Opposite to the side of the 
triangle whose length depends on, in addition to the CKM parameters \V c b\ and \V u d\, the 
elements involving heavy-light quark mixing \V u b\ and \V C d\, the angle f3 is presently well- 
measured. So precision determination of them is central to the unitarity testing. Exclusive 
processes offer an indispensable avenue to understand these parameters. The decays of heavy 
mesons into a light pseudoscalar meson plus an electron and its antineutrino can proceed 
at electro-weak tree level and are much less sensitive to new physics, and accordingly they 
could serve as preferred exclusive channels to probe both elements that we take interest in, 
namely, \V u b\ and \V c d\- Then we are confronted with calculation of the hadronic matrix 
elements, say, that for the B° — > 7r + transition parameterized usually as 

(7r(p)|« 7M 6|S(p + q)) = 2f^(q 2 )p, + (f^(q 2 ) + f B ^(q 2 ))q„ (1) 

with the momentum assignment specified in brackets, and f+^ n (q 2 ) and f^^iq 2 ) being 
the form factors describing QCD dynamics in the decay, of which only the former is related 
if the small electron mass is neglected. Combining the partial rates measured in some q 2 
bins with the form factor predictions of different QCD approaches, one could achieve the 
values for related \Vij\. Another approach is to fit the experimental observations using 
the various form factor parameterizations. In such way, a strong constraint is imposed on 
q 2 distributions of the form factors such that one may obtain a precise estimate of the 
products / + (0)|Vy |, in which case theoretical task boils down to estimating the form factors 
/ + (0) at q 2 = 0. Requiring a good knowledge of the form factors, the exclusive avenues 
to \Vij\ are theoretically more challenging than inclusive approaches. The continual data 
updates have aroused one's enthusiasm for exploring heavy-to-light transitions to approach 
an understanding of the CKM parameters. In the wake of the recent accurate measurements 
of the semileptonic processes by the BaBar [l|, 0] and CLEO [3, 0] collaborations, new 
progress has been achieved in this respect. Some extent of tension, however, still holds 
between inclusive and exclusive extractions of \V u b\. A global data- fitting from CKMFitter 
5| and UTfit (J is in favor of a smaller \V u b\ than inclusive determinations. One can 
be referred to 7fl for a comprehensive overview of the current status of the CKM matrix 
elements. 

Developed from QCD sum rule technique, light cone sum rules (LCSRs)^, §] have be- 
come a powerful competitor in making predictions for heavy to light transitions. Comple- 
mentary to lattice QCD (LQCD) simulations, this approach is successfully applied to study 
B decays 9-[l6]: whereas the former are available for the high q 2 , LCSR calculation is appli- 
cable for the low and intermediate q 2 . Utilizing the LCSR predictions for f^ n {q 2 ), one has 
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launched a painstaking investigation into \V u b\ [lll-ll4|. with a consistent result with those 



using LQCD. The same approach has also been taken to understand D — >■ tc,K decays in 
0,0, the resulting sum rules [l8| being employed to extract \V cc i\ and \V CS \ . 

The uncertainties in the light meson distribution amplitudes (DAs) involved in the sum 
rules, however, would have different degrees of impacts on the results. To gain enlightenment 
on how to further improve accuracy of the LCSR calculations, it is essential to look into the 
role played by each of the higher twist DAs. A systematic numerical analysis shows that 
whereas the twist-4 effects account for only a few percent of the total sum rule results, the 
chirally enhanced twist-3 contributions are numerically large enough to be comparable with 
the twist-2 ones in the B meson cases, and even about twice as large as the latter for D 
decays. As a result, there are a few problems left unsolved. To start with, one might doubt 
whether the potential twist-5 effects are negligible in particular while assessing D decays. 
Secondly, the sum rule pollution by twist-3 would be serious on account of the combined 
uncertainties of the DAs and chiral enhancement factor. Finally, since there is an extremely 
different sensitivity to twist-2 between the sum rules for /^"^(O) and f+~*"*(0), a successful 
LCSR application to the latter does not necessary assure, with the same inputs, a reliable 
LCSR prediction for the former. For the moment, these issues are difficult to essentially 
settle within the LCSR framework. The trick suggested in (9), [l5| is available as a temporary 
scenario to approach them. 

Focusing on f^ n (0) and f+~*"*(0), in this work we intend to reconsider heavy to light 
transitions in the revised LCSR version so as to provide a calculation independent of the 
traditional LCSR ones and further a determination of the associated CKM parameters. We 
will expound that this approach does not involve the twist-3 and-5 DAs, which enables us 
to get an understanding of the form factors to twist-5 precision only resorting to the known 
twist-2 and -4 DAs and to perform a cross-check between the resulting LCSRs for f+^"*(0) 
and /+~ +7r (0). This paper is organized as follows. In the following Section we put forward 
our derivation of the sum rules in question, including a detailed next-to-leading order (NLO) 
QCD calculation in twist-2 approximation, and elaborate on the key technical points. The 
modifications and improvements made are also addressed in comparison with the previous 



calculations [15l. |16|. In Section 3, after discussing assignment of the parameters for which 
updated and consistent findings are selected as inputs, we shift into numerical computation 
with a systematic error discussion included, by means of up-to-date experimental data, and 
present our LCSR results for f+~* n (0) and /+ ^""(O) and the determination of \V u b\ and \V c d\- 
Too we report on an estimate of the decay constant /d, as a by-product. The final Section 
is devoted to a summary. 
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II. QCD CALCULATION OF f^(0) AND f^ n (0) 

The starting point of LCSR calculation is to consider a correlation function with T 
product of currents sandwiched between the vacuum and a light meson state L. In the 
coordinate space and for large and negative virtuality of the current operators, the correlation 
function can be in form expanded, in the small light cone distance x 2 ~ 0, as, 

correlation function ~ C m (x)(L(p)\O m (x, 0) 1 0) , (2) 

m 

where C m (x) are the Wilson coefficients, O m (x, 0) the nonlocal operators built out of quark 
and/or gluon fields, and the matrix elements {L(p)\O m (x, 0) 1 0) have an expansion form in 
term of the light cone DAs \E f( - n ^ with increasing twist n. The power series *YlC n (x ■ p) n 
(x ■ p ~ 1 for a large external momentum p) appearing in the expansion process are summed 
up effectively, which works out some of the problems with the expansion in the small distance 
x ~ 0. Switching (2) to momentum space, we have 

correlation function ~ V^^ir^ <S> (3) 

n 

a factorized form with the hard kernel being convoluted with vl/W. Whereas the process- 
independent \I/( n ) parameterize the long distance effects below a factorization scale /i, the 
process-dependent amplitudes describe the hard-scattering dynamics above /i, which 
are perturbatively calculable and have the following expansions in a s : 

Ti n) = T (n) + ^F T (n) + . . (4) 

If calculation is restricted to 0(a s ) accuracy, we need just to estimate the leading order (LO) 
contributions and NLO corrections . Then the remaining procedure is standard. 

Now let us take up our LCSR calculations of f+~* n (0) and f^ n (0). Allowing for the 
similarity of the two situations, for definiteness we would like to concentrate on the former. 
Moreover, throughout the paper the chiral limit m w = is taken. We follow 0, [3] and 
adopt the following correlation function: 

n„(p,g) = i f d 4 xe^(n(p)\T{jV+ A (x),jr +s (0)}\0) 

= F((p + q) 2 )p, + F((p + q) 2 )q,- (5) 

Here we substitute the chiral currents J^ +A (x) = ■u(x)7 Ai (l + 7 5 )6(x) and J^ +s = 
rribb(0)i(l + 7s)<i(0), respectively, for the operators adopted usually J^(x) = u(x)^ '^(x) 
and Jb = mJ}{fS)i^d{Q) . The operator replacements do not violate renormalization group 
invariance of the correlation function, for both J^ +A and j£ , like the latter two, have an 
anomalous dimension of zero, and however make the correlation function receive an addi- 
tional contribution from the set of scalar B mesons. In view of that the mass of the lowest 
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FIG. 1: Tree-level Feynman diagrams contributing to the correlation function. 



scalar B meson is slightly below the one of the first excited state of the pseudoscalar B 
mesons, we could safely isolate the pole term of the pseudoscalar ground state from the 
contributions of higher resonances and continuum states. 

For the present purpose, it is sufficient to consider the part proportional to p^ in (5), that 
is, the invariant function F((p + q) 2 )- It has the pole term of interest to us, 

Fpoieiip + q = — 2 1 — , — vT' 6 

m B - [P + Q) 

where rns and fs indicate, respectively, the B meson mass and decay constant defined as 

(S|W 76 d|0> = (7) 

The spectral function p H {s) is introduced to reckon in the higher state contributions in 
a dispersion integral starting with the threshold sjf, which should be assigned near the 
squared mass of the lowest scalar B meson. At this point, what remains to be done is 
the light cone expansion calculation on F((p + q) 2 ), from which the corresponding QCD 
spectral function (s) is extracted in order to get the sum rule for f^ n (0) by matching 
the Borel improved theoretical and phenomenological forms with the duality assumption 
p H(s) = p QCD (s)0(s- 3%). 

The light cone expansion of (5) goes effectively in the large space-like momentum region 
(p + q) 2 — ml « for the bd channel. At tree-level and to NLO in the light-cone expansion 
of the b quark propagator, it can be illustrated by the two Feynman diagrams as depicted 
in Fig.l. In comparison with Fig. 1(a), which corresponds to the leading term in the quark 
propagator and illustrates the two-particle contribution, Fig. 1(b) portrays the three-particle 
Fock state effect due to the soft-emission correction to the free quark propagator, which is 
expressed as 



(8) 
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The contribution of Fig. 1(a) to F((p + q) 2 ) is easy to estimate, using the definition of the 
pionic two-particle DAs: 

U f 1 



(n(p)\u a (x)clp(0)\0) x ^ o = i^J due^ x [(jh,)^{u) 



1 1 r u -I 

+ J^(Jl5)f3 a X 2 <fr<l7r(u) -i-(^ 5 )l3a J ^{v)dv , (9) 

where u is the fraction of the light cone momentum p + p 3 of the pion carried by the 
constituent u quark. While <f n (u) denotes the twist-2 DA, both <p^ n (u) and <p^ n (u), which 
are accompanied by the chiral enhancement factor fi n , have twist-3, and the other two 
functions are both of twist-4. From the following trace form, which emerges obviously as 
one works in the momentum space, 

Tr{[d(up)u(up)} 7m (1 + 75) {d + up 1 + m b )(l + 75)}, (10) 

N v ' 

wavef unction 

we see readily that the twist-3 components make a vanishing contribution to the light cone 
expansion, because of the corresponding Dirac wavef unctions. In fact, the same happens to 
the three-particle situation, as shown from a straightforward computation with (8) and the 
decomposition: 

(7r(p)\u a (x)g s G^(vx)dp(0)\0) x ^ = l - j Va^ x ^ +a ^[ifU° pX ls)p a 

+ P P (PnX„- Pu X u ) (<M ^ + j _ ik €lwsMfia 

x { (p x g 5p - p 5 g xp )*M + pP{p5x *- pXx5) (^(o.) + y^m)) }] , (11) 

where is the gluonic field strength tensor and T>a.i = daida2da^5(l — a x — a 2 — ^3); 
$ 37r (o;j) indicates the twist-3 component of the three-particle DAs, and the remaining func- 
tions are all of twist-4. In the usual LCSR application to D decays, the chirally enhanced 
twist-3 terms provide a leading contribution, which engenders much negative influence as 
af orement ioned . 

At present, the two-particle contribution F^ 2p \(p + q) 2 ) can be written down in a form 
that the DAs are convoluted with the corresponding LO hard scattering amplitudes, 

F Q {2p \{p + qf) = ~U fdu [T^\{p + q) 2 ,u)vM 

Jo 

- T W ((p + q) 2 ,u) / Mv)dv-T^\(p + q) 2 ,u)Mu)l (12) 
Jo 
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with 



T {2 \(p + q) 2 ,u) = -2 



4 A \{p + q)\u) = 2- 



u 

u 2 d 2 



d 

u— + 1 



2(ml — u(p + q) 2 ) du 2 
The three-particle contribution is of the following convolution 



(13) 
(14) 
(15) 



F {3p \(p + qf) = -U fdu T { *\(p + q) 2 ,u)lUu) 

Jo 



with 



"(4) 



t-((p+?)» = 2- 



u 



d 



ml — u{p + q) 2 ) du 



(16) 

(17) 
(18) 



and 



hJu) 



da i 



<£V(«i) 



dv 

(«-ai)/(l-oi) U 

(a 



oi3=(u— ai)/v 



(19) 



Then we can attain the imaginary part of Fq CD ((p + q) 2 ) = F^ 2p \(p + q) 2 ) + F^ 3p \(p + q) 2 ) 
via estimating the ones of the hard kernels in (13-15) and (17), and further the desired QCD 
spectral function Pq CD {s). The result is as follows: 



P Q » CD {s) 



f 1 s 
2/,r / duS (l — u — 

Jo m b 



u 





1 d 2 



u ( d 

iPirW + — \U— + 1 



4ml du 2 



b 

4n(u) \—I^{u) 

m{ du 



mi 



du 



1p4n(v)dv 



(20) 



At twist-4 level, we have provided a complete LO light cone QCD representation for 
F((p + q) 2 ). In its present form, the ensuing continuum substraction could be enforced 
systematically for the twist-4 as well as twist- 2 parts, with the known QCD spectral function. 
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16j in which the twist-4 contribution is 



This improves explicitly the previous treatment 
written down in a form not suitable for continuum substraction. 

Our main task is to evaluate the gluon emission effect on the LCSR for f+^"*(0) at one 
loop level. It should be sufficient for this purpose to calculate the NLO parts of the leading 
twist-2 and the chirally enhanced twist-3 contributions. To be specific, we are about to 
compute the six Feynman diagrams plotted in Fig. 2, to the accuracy in question. Fig. 2(a) 
depicts diagrammatically the hard-exchange correction between the outgoing and spectator 
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FIG. 2: One-loop Feynman diagrams contributing to the correction function. 



quarks in the B — > tt transition. From the nature of the correlation function, we deduce 
easily that there is no UV divergence in Fig. 2(a) or it could not be canceled out. Of the other 
figures, Figs.2(b, e) and Figs.2(c, f) involve, respectively, the partial one-loop contributions 
to the I + 75 vertex and to the 7 M (1 + 75) one, while Fig.2(d) does the remaining loop 
contribution to both operators. It is conceivable that each of these five includes both UV 
and IR divergences, except Fig. 2(d) which is merely UV divergent because obviously if any 
IR divergence arises it can not be reasonably absorbed into a pionic DA. 

It is found that the twist-3 components still produce no effect at one-loop level, for the 
same reason as in the tree-level case. Hence the NLO computation is reduced to a calculation 

of the 0(a s ) correction to the LO twist-2 contribution (for brevity, hereafter we indicate 

(2) 

the LO twist-2 contribution by the symbol T instead of T and the corresponding NLO 
correction by T 1; up to a prefactor a s /47i). We work in the Feynman gauge. In addition, 
we use the dimensional regularization and MS scheme to deal with the ultraviolet (UV) 
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and infrared (IR) divergences appearing in the calculation, such that the LO evolution 
kernel of [19j achieved early in the same prescription is available for a proof of QCD 
factorization for the resulting twist-2 contribution to F((p + q) 2 ) as we attempt to segregate 
the long distance contribution from the perturbative kernel. The calculation is tedious and 
complicated. Here we present, for the first time, some details of the diagram calculation. We 
summarize the divergence contribution to T\ from each of the diagrams in Fig. 2 as follows, 



ur' 



r)ln(l — r) 



r ) ln(l — ur) 



Am, 



1 — ur 
1 

1 — ur 

2 + ur 
(1 — ur) 
2 



ur 



ur 



+ 1 1 A 



^IR 



2A 



uv 



ur 



2 Auv, 



-ln(l -ur) + l) A IR - -A 



T$ +f) (u,r) = - 

K " I — ur 

Here u = 1 — u, r = (p + q) 2 jm\ and 



A 



IR 



-A 



ur 



uv- 



Air (A 



uvj 



1 



e lR \ £ VV 



- 7£ + ln47r 



(21) 
(22) 
(23) 

(24) 
(25) 

(26) 



with the Eyy and £ir introduced to regularize the UV and IR divergences, respectively. 
Obviously, the yielded results expected. 

Adding all the divergent and finite terms together, we have the NLO correction 

1 — ur\ 1 — r — uur 2 



Ti(u,r) = 2 
+ 



1 — ur 
Qur 



1 / 3-2 1n(l-r) 1 - / ' ~ " r 



2 In 



+ 2 

+ 2 
+ 

- 2 



1 — ur 
1 



1 + ur 

-A uv + — 

z 1 — urY 



1 



1 



ur r(l — ur) \ur 2 
1 1 



3 - 3 ln^ + 1 
u z ur 

u 2 



1 



ur ur 2 



ur 
2 



ur + u r + u 
H ~, 2 



1 — ur 
(ln 2 (l-r) + Li 2 (r)) 

2 1 — ur 



1 — r 



ln(l — r) 



uur' 



A 



IR 



1 — ur 



u[ur ) 
L — ur 



uur" 



1 — ur 



uur z 



h4 



ln(l — ur) 



(ln 2 (l — ur) + Li 2 (ur)) 



(27) 



ln(l-f) 



with the dilogarithm Li 2 (x) = — dt 

Keep in mind that up to now the quark mass has been treated as a bare quantity. A mass 
renormalization must be performed in the MS scheme, in order to have a UV renormalized 
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hard-scattering amplitude T via adding T\ to To. It can be done by making the parameter 
replacement m& — > Z m m b m the related expressions, with the renormalization constant 
Z m = 1 — 3 A uv a ^ F ■ As a result, the tree level expression (13), to the accuracy required, is 
modified to the form 



T (u,r) 



a s Cp Ylur 



uv, 



(28) 



ur — 1 47r (1 — ur) 
but the NLO term T\(u,r) keeps its form unchanged. Here m& entering r should be under- 



stood as the MS mass. The additional UV divergent contribution in (28), as it should be, 
precisely cancels out the one of (27). Then a complete UV renormalized result is obtained 

as 

OL„Cp. 



T(u, r) 



T (u,r) + 



An 



-Ti(u,r) 



+ < 

+ 2 
+ 2 
+ 



1 — ur 
1 + ur 
1 — ur) 2 
1 



3-2 ln(l -r 



1 — r — ur 



2 In 



ur^ 



1 — ur \ 1 — r — uur z 



1 — r 



uur z 



A 



IK 



3-31n4 + J-) 

fj, z ur J 



ur r(l — ur) 

1 1 
1 — ur ~ n " 2 

+ 



1 



In 



ur' 



4 ur + u 2 r + it 



Mr 2 ]_ _ ur 
(ln 2 (l-r) + Li 2 (r)) 

2 1 — ur 



ln(l — r) 



1 — wr 
2 



m(ot) 2 
1 — ur 
uur 2 



- 2 



uur 2 



1 — ur 

(ln 2 (l - ur) + Li 2 (wr)) }>. 



In^' 



ln(l — ur) 



1 — ur 

We need to add that superior to use of the pole mass for the b quark 



(29) 
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employing 

the MS mass could render not only the calculation free from some element of uncertainty 
but the physical meaning more obvious even when the calculation is performed at QCD tree 
level, as shown in (28). 

To proceed, we embark on handling the IR divergence term, 



T IR (u,r) 



2A 



IK 



ur l 



3 - 21n(l - r) 



' ~ Ur - - 21n(l 



ur , 1 — r — uur 2 



ur' 



r 



uur 2 



If we try to subtract the divergent part from the UV renormalized hard amplitude to repre- 
sent the invariant function F((p + g) 2 ) in the form of QCD factorization, it has to abide by 
the form 

• i 



T lR (u,r) = - A 



IK 



dv V (v,u) T (v,r), 



(30) 



where Vq(v, u) is the kernel of the evolution equation of the pionic twist-2 DA 19J|. As checked 
readily, this is indeed the case. We can therefore eliminate the divergence by defining a scale 
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dependent DA as 



ip*{u,n) = (p v (u) - Air 



An 



f 

Jo 



dv V (u,v) (p w (v), 



(31) 



which is convoluted with the perturbative kernel T H (u,r, /i) = T(u,r)\A m =o- As a result, 
the twist-2 contribution to F((p + q) 2 ) observes, at NLO, the following QCD factorization: 



F QCD ((p + q) 2 ) = -U [ duT H {u,r^) <p v {u,fj,). 

Jo 



(32) 



Up to higher order corrections in a s , /x dependence of T H (u,r,fj) compensates that of 
ip w (u,[/). It should be understood that in the above operations the factorization and renor- 
malization scales have been set identical for simplicity. 

Having in hand the hard kernel available, we can calculate the QCD spectral function to 
write F QCD ((p + q) 2 ) as a dispersion integral. For r = (p + q) 2 /m 2 = s/m 2 > 1, we have 



1 r 

p QCD (s) = f w / dr] lmT H (u,r,/j,) <p v (u, n)\ u =r,/r, 



(33) 



2vr 



lmT H (u, r, jx) 



u=r//r 



6 - 31n^ 
A* 



7 -ti 2 - 2Li 2 (l - r) + 21n 2 (r - 1) - 2 Mnr + - - 1 ) lu(/ I ) 



- 2 4 - 31n 



A 4 



1 + 



drj 



21n(r - 1) 1 - In 



m 



+ 2(9(77 - 1) 



, mt 
-2 +ln^ 

fi 2 

1 



"4 ln(?7 - 1) 


1 

+ 7 


i] — 1 


rj-1 



In 



+ 



(r — l) 2 r 



r(r — rj) 
+ 29(1-7]) 

1 



1 — r 

7]T 



lnw — 2 In 



In 



mi 



(rj-1)'' 



rj — 1 
V - 1 



+ 1 - In— f- 
fi 2 

Inrj i](r — 1) — r 



77-1 

,2 



2r?7 2 



r 1 mA 1 
In, ^ + - - ln-£ 



(r — l) 2 / 



+ 



ln- 



r(r — rj) \ (r — 1) 



/i 2 



1 -r 

r(r — 77) 



(34) 



where we take the operation 



F( V ) 



1-77 



+ 



F(y)-F(l) 
I-77 



(35) 



to avert the redundant divergences possibly occurring as the integral in (33) is performed 
over the interval [0, r]. 
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Using (33) and counting the twist-4 contribution covered in (20), we have the final sum 
rule for the product /b 1+^(0) 

/ B /r-(o )e -* = -fh I**.*-* 1 - /• / "' l *,i m r(^„,4,,) ^(<^) 

2nm B J m 2 s J V s m 2 b ) s 



U f 1 , -J*L ( ud 2 (p i7r (u) 
+ — w \ due 



+ uip ilT {u) + / dvif) in (v) - ^^{u] 
= K(s^M 2 ), (36) 

with M 2 indicating the Borel parameter with respect to (p + q) 2 and Uq = m 2 /sQ. 

Converting (36) into the corresponding sum rule for D — > n transition by a simple re- 
placement of the parameters, we put an end to our derivation of the LCSRs for f+^"*(0) 
and /^^(O), to 0(a s ) precision in twist-2 approximation and at tree-level for twist-4 con- 
tributions. 

We close this Section with a few remarks. Albeit the LCSR calculations are done at 
twist-4 level, the results remain valid to twist-5 accuracy. The reason is simple. The twist-5 
DAs as well as twist-3 ones play no role in the present context due to the Dirac structures of 
the related nonlocal operators, of which both d^^u^O) and d(x) 07^75^(0), as sandwiched 
between the vacuum and a pion state, bring about a chirally enhanced twist expansion. 
Apart from helping reduce sum rule pollution by long-distance parameters, the disappear- 
ance of twist-3 and -5 components from the light cone expansions guarantees the resulting 
LCSRs well convergent. We are going to return to this point in the following Section. 



III. CHOICE OF THE INPUTS AND NUMERICAL DISCUSSION 

Presently, theoretical estimates of f+^(Q) and f+^(0) with twist-5 accuracy are ob- 
tainable in the sum rules to have been given and the inputs to properly be selected. On the 
experimental side, from the measured shapes of the form factors for B — > tiIu, the CKM 
matrix element \V u b\ multiplied by f+^ n (0) is numerically inferred as[l|: 

f+^(0)\V ub \ = (9.4 ± 0.3 ± 0.3) x 10" 4 . (37) 

For the semileptonic processes D — > 71W, a similar manipulation |3] gives 

7+^(0)1^1 = 0.150 ± 0.004 ± 0.001. (38) 

Then the yielded theoretical predictions could have \V u b\ and \V c d\ extracted from these 
up-to-date data. 

Aimed at determining \V u b\ and \V c d\, we must do our best to enhance reliability of the 
LCSR assessments for the form factors in question. So special care should be taken when 
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making our choice of the parameters entering the sum rules. The main sources of uncertainty 
are, of course, the related DAs, which can merely be understood at a phenomenological 
level. Based on the conformal symmetry of massless QCD, we can parameterize these DAs 
by expanding them in terms of matrix elements of conformal operators. The twist-2 DA 
(Pn(u) is of the following expansion in the Gegenbauer polynomials: 



6uu ^1 



a2(yu)Cf^ 2 (u — u) + a±{n)C^/ 2 (u — u) + ■ 



(39) 



with the even moments a 2ri (/z) remaining to be determined. The Gegenbauer polynomials 
of higher-degree (large n) are rapidly oscillating and so one neglects usually their effects on 
the numerical integrals included in the sum rules by retaining only the first few terms of the 
expansion. Some scenarios have been put forward to examine the higher-moment effects. 
We are willing to mention the prescriptions suggested in j2o| and in 16 1. In 20] Ball and 
Talbot (BT) presume that a 2n fall off as powers of n, a 2n oc l/(n + l) p , in order to build 
a DA model. In comparison, authors of 16| consider a modified transverse momentum Ki 
dependent Brodsky-Huang-Lepage (BHL) wavefunction, 



V v {u,K±) = [l + B n C2 /2 (2u-l) + C n Cl /2 (2u~l)] 



At, 



ml — u 



-exp 



K 2 



m. 



ml — u) 



(40) 



which is integrated over |Kj_| < \i to give a twist-2 DA. Phenomenological studies with 
both models are in support of the rationality of using an expansion truncated after n = 2. 
We stick to such disposal. In one-loop approximation taken as default for all the renormal- 
ized parameters except QCD coupling, a 2 (/x) and a 4 (/i) respect the renormalization group 
equations 



25C F 



a 2 (/z 2 ) = [L(fi 2 , fi>i)] 6/3 ° a 2 (/ii), 



(41) 



91CV 



04(^2) = [L(// 2 ,/xi)] 15,30 a 4 (/ii 



(42) 



with L(/i 2 ,/ii) = asfa) ' = 4/3 and f3o — 11 — ^j-, n/ being the number of active quark 
flavors. To our knowledge, all the existing estimates for a 2 (/i) are basically consistent with 
each other and have the averaged central value of 0.25 at fi — 1 GeV. In the light of 
the current experimental constraints imposed on LCSR calculations, a 2 (l GeV) appears to 



prefer varying between 0.16 — 0.19 [12H14J. l2ll |. The situation is not optimistic about a^fi) 



The findings differ among the various studies to a large extent, and even there would be 
a difference in sign between numerical estimates. Fortunately, the sum rule results depend 
less sensitively on 04 (/i) than on a 2 (/i). We would like to use as a consistent input the 
findings 14], a 2 (/i = lGeV) = 0.17 ± 0.08 and a 4 (/i = lGeV) = 0.06 ± 0.1, from fitting 
the LCSR calculation of the pionic electromagnetic form factor to the recent experimental 
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observation. Concerning the twist-4 DAs, the three-particle components are specified by 
only two parameters to NLO in conformal spin, and are of the following forms: 



$47r(aj) = 1205le n (a>i - a 2 )aia 2 a 3: 



(43) 



■ — I205laia 2 a 3 
*4x(ai) = 30^a 3 2 (l - a 3 ] 



$4n(&: 



3+2^(1 



e n (l - 3a 3 ) 



2a, 



+ 2eJl - 2a 3 



(44) 
(45) 
(46) 



where the nonperturbative quantities 5% and e n have the scale dependence 



SCp 

5l(fi 2 ) = [L(n 2 ,n x )Y^5l(iix), 
(%e«)M = [L(m 2 ,Mi)]%X)(mi), (47) 

and the parameter values [22f 6% = (0.18 ± 0.06)GeV 2 and e n = fu 47T (u 47T = 0.2 ± 0.1) 
normalized at 1 GeV, which are to be adopted as inputs. Resorting to equation of motion 
the two-particle components, without introducing any new parameter, can be understood 

as 



^A-k(u) 



= -^-5lu 2 u 2 + 86le n {uu(2 + 13mm) 

+ 2m 3 (10 - 15u + 6u 2 )\nu + 2m 3 (10 - 15u + 6u 2 )\nu}, 



(48) 
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The MS quark masses m b and m c comp 



(49) 

y with the proverbial LO evolution equations, 
sum rule results with four-loop precision, 
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The bottomonium 23( and charmonium [23 
m b (rh h ) = 4.164 ± 0.025 GeV and rh c (m c ) = 1.29 ± 0.03 GeV, are applicable well to the 
present discussion. As far as QCD coupling goes, we use two-loop running down from 
a s (M z ) = 0.1176 ± 0.002 [251 ] . Additionally, the factorization scales are assigned, according 
to the typical virtuality of the heavy quarks, as fib — 3 GeV and u c = 1.5 GeV in the 
respective cases of B and D mesons. 

Among the hadronic parameters are the decay constants fo, and /„., apart from 
the heavy meson masses determined experimentally j25| as tub = 5.279 GeV and rri£> = 
1.865 GeV. The value of f n is measured at f n = 130.4 MeV [25j, from the exclusive 



processes it — > fiu^ and tt — > fw^l- Recently, an updated measurement of fo has already 
been reported by the CLEO collaboration [4], f& = 205.8 ± 8.9 MeV. However, it is on the 
basis of combining the experimental data on fo multiplied by \V c d\, 



f D \v c 



cd 



46.4 ± 2.0 MeV 



(50) 
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and the assumption \V c d\ = \V US \ = 0.2255±0.0019, and hence could only serve as an input in 
the sum rule calculation of /+ "^(O). Instead of a direct estimate of /+ "^(O), we consider the 
sum rule for the product /dZ+^CO), which in conjunction with the experimental numbers 
(38) and (50) allows us to consistently make predictions for the quantities f+~* n (0), \V c d\ and 
fo as well. By contrast, leptonic B decays are made difficult to detect experimentally by 
higher helicity suppression. To have a measurement analogous to (50), the only opportunity 
is furnished by B — > rv T well established lately j26|. Nevertheless the results yielded in the 
SM are less persuasive. The reason is that these modes turn out to be sensitive to possible 
extensions of the SM such as the two-Higgs doublet models and minimal supersymmetric 
extensions. We must have recourse to theoretical predictions for to make an assessment 
of /^^(O). As a consistent choice, here we make use of the interval Jb = 214+^ MeV [3] 
from a sum rule with the MS quark mass. 

The remaining parameters are intrinsic to the sum rules, containing the effective threshold 
Sq (sq) and Borel variables M 2 . The former can be set at the neighborhood of the squared 
mass of the lowest scalar B meson (D meson). An alternative manner, which has proven to 
be more effective, is through use of an auxiliary sum rule obtained, for example, by taking 
logarithmic derivative of 1/M 2 for (36), 

m| = -^lni^,M 2 ). (51) 

Requiring the measured value of the B meson mass to be reproduced precisely from the above 
sum rule, we get the effective interval = (34 ± 0.5)GeV 2 in accordance with the sum rule 



estimate in heavy quark effective theory 271 ] . Similarly, is fixed at (6.5±0.25)GeV 2 . The 
Borel intervals could be specified in the standard procedure. We have M 2 = (18 ± 3) GeV 2 
and M 2 = (6 ± 3) GeV 2 , corresponding to, respectively, the sum rules for B and D mesons. 
As both inherent parameters vary within their separate ranges allowed, it is demonstrated 
that the twist-4 effects are kept at a numerical level less than 4%, and also the continuum 
contributions are highly suppressed, not exceeding 20%. 

Using the inputs given above, the numerical discussion can be done. Our sum rule result 
for f B f+^(0) reads 

/s/+^(0) = 59±J° MeV, (52) 

with the uncertainty achieved by adding in quadrature all the errors caused by variations 
of the inputs, of which the scale parameter \i is set to the interval between (2.5 — 6.0) GeV. 
We address this result is because it is independent of the value for fs and therefore of less 
uncertainty, and moreover is convenient for a numerical update of the sum rule for f?^ n (0) 
once the theoretical estimate of fs gets improved in the future. Substituting the parameter 
value for into (52), we obtain 

/f-*(0) = 0.28t° : °2- (53) 
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FIG. 3: Dependence of the LCSR for fg n {0) on the Borel parameter M 2 (a) and on the factorization 
scale fi (b) 



TABLE I: The LCSR result for f+~*"*(0) with the uncertainty estimates due to the variation of the 
input. 



Central value 


M 2 


s- 8 
s o 


1* 


m b 


fs 


a\ 


a\ 


/f-^(0) 
0.277 


+0.002 


+0.007 


+0.05 


+0.008 


+0.007 


+0.008 


+0.01 


-0.001 


-0.008 


-0.01 


-0.008 


-0.009 


-0.008 


-0.01 



Illustrating stability of the numerical result, we display the variations of the sum rule for 
f?~**(Q) with the Borel and the scale parameters, respectively, in Figs. 3(a) and 3(b). It 
is distinctly observed that the M 2 dependence is considerably weak in the Borel interval 
required, and there is a moderate \i dependence. Furthermore, to have an explicit under- 
standing of the role that every source of uncertainty plays in the uncertainty evaluation, 
we collect in Tab.l the individual uncertainty contributions estimated by altering each of 
the inputs within its specified range. Those not listed therein are tiny and included in 
the total uncertainty. A comparison is drawn among the LCSR predictions for ff~*"*(0) in 
Tab. 2, there being a result quite close to one another. We can understand it as follows: 
(1) No matter which of the two correlation functions one adopts for a LCSR estimate of 
that quantity, the light-cone expansion reveals a good convergence, as will be addressed. (2) 
All these calculations employ essentially the same inputs for the leading twist-2 DA, along 
with a fg consistently determined from the sum rules. It is exceptionally hard to have a 
LQCD calculation to compare with, since the pionic energy goes beyond the restriction by 
the lattice spacing. Nonetheless, it is claimed [28] that f+~^ n (Q) is estimable in an improved 
LQCD simulation, with the result /£ ^(0) = 0.27 ± 0.07 ± 0.05. 

Now the experimental measurement (37), with the aid of the theoretical prediction (53), 
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TABLE II: Comparison of theoretical predictions for the form factors f+^ n (0) and /_P~ s>7r (0). 



Approach [Ref.] 


/f-^(o) 


/f^(0) 


LCSR [14] 

\m 

[lZj 
[13] 
[18] 
This work 


0.281 ± 0.05 
0.258 ±0.331 

26+ - 04 
u - ZD -0.03 

n 9 o+0.05 


0.63 ±0.11 

u -°' -0.07 

0.62 ± 0.03 


Lattice QCD [30] 
[29j 
[28] 
[31] 
[32] 




0.57 ±0.06 ±0.02 
0.64 ±0.03 ±0.06 
0.74 ± 0.06 ± 0.04 

0.666 ±0.029 
0.65 ±0.06 ±0.06 



allows for extracting the desired CKM matrix element \V u b\. We have the interval: 

\V ub \ = (3At 2 6 ± 0.1 ± 0.1) x 10- 3 , (54) 

where the first error originates from the uncertainty of f+~* n (0) and the others do from 
the corresponding experimental ones. Obviously, an analogous result can be extracted in 
the other LCSR estimates of /^^(O) in Tab. 2. There is also a similar determination from 
matching the LCSR calculations and experimental partial rates for q 2 < 12GeV 2 [14]. All 
these are upheld bv the findings in LQCD simulations for a high q 2 and consistent with the 
CKM fit upshots 

Corresponding to (52), the product /d 7+^(0) has the numerical value 

fof+^ = H7l? MeV. (55) 

Tab. 3 provides a summary of the major uncertainty contributions to the sum rule. As 
exhibited in Figs. 4(a) and 4(b), the stability of the sum rule holds as well as in the B meson 
situation, as M 2 changes in the interval specified and /i ranges from 1 to 3 GeV. Intriguingly, 
using the same inputs as ours for most of the parameters this quantity is explored in the 
LCSR approach [l8[ and the yielded result /^/^(O) = 1371^ MeV is compatible with our 
prediction within the errors, but showing a larger central value. We remark on this difference. 
The twist expansion in x 2 is the basic thought of the LCSR approach. For heavy to 
light transition, such an expansion must match the one in the inverse of heavy quark mass 
rriQ. One shows, indeed, that in the heavy quark expansion the end point behaviors of the 
higher-twist DAs entering a traditional LCSR might modify, but does not violate the twist 
hierarchy. For instance, the twist-3 term, which is formally 1/toq suppressed versus the 
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FIG. 4: Dependence of the LCSR for /d/o~(0) on the Borel parameter M 2 (a) and on the factor- 
ization scale \i (b). 



twist-2 part, behaves the same as the latter in the heavy quark limit. However, an explicit 
calculation with a finite mq demonstrates that whereas the twist expansion works better 
for B decays, there is a considerable numerical violation of the hierarchy relation in the D 
meson cases, where the twist-3 components contribute to the sum rules much more than the 
twist-2 ones due to the chiral enhancement factor fi n > 1. The fact that the sum rule for 
f+^iO) is poorly convergent implies that the twist-5 effect is not negligible and should be 
considered, even if we work in twist-4 approximation. Currently nothing is known, however, 
about the twist-5 DAs except that they provide the sum rule with a term formally l/mj? 
suppressed with respective to the twist-3 one. To have a sketchy understanding of their 
influence on the LCSR calculation, authors of [18| suppose that the ratio of the twist-5 to 
-3 parts is identical to the one of the twist-4 and -2 terms, while in 17| the twist-4 term 
is multiplied by a factor of 3. Anyway, it is still obscure that how much on earth do the 
twist-5 components, in particular those with the chiral enhancement factor, contribute to 
the sum rule for f+^"*(0). We leave it as an open question until a reliable twist-5 model 
wavefunction is presented. Given that the present scenario ensures, to twist-5 precision, the 
light-cone expansion to converge well whether for B or D decays, this issue gets, at any rate, 
settled provisionally. 

Let us go back to our numerical calculation. Combining the sum rule prediction (55) 
with the product of the two experimental numbers (38) and (50), we could yield the square 
of \Vcd\ and further the magnitude of V c d- 



IK 



cd 



0.244 ± 0.005 ± 0.003 ± 0.008, 



(56) 



where the first and second errors are of an experimental origin and the third is due to the 
theoretical uncertainty. This result deviates by about 2% from the Wolfenstein approxima- 
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TABLE III: The LCSR result for fDf+~*™(fy with the uncertainty estimates due to the variation 
of the input. 



Central value 


M 2 


s 





m c 


a\ 


a\ 


U3\ 


51 


0.117 


+0.0007 


+0.0015 


+0.0062 


+0.0001 


+0.005 


+0.0017 


+0.0003 


+0.0014 


-0.0001 


-0.0016 


-0.0049 


-0.0003 


-0.005 


-0.0016 


-0.0002 


-0.0013 



tion \V cd \ = \V US \= 0.2255 + 0.0024 [25], and is in good keeping with \ V cd \ = 0.234 + 0.007 + 
(1002 ± 0.025 extracted from (38) by using the LQCD estimate f£ n (0) = 0.64 ± 0.03 ± 0.06 



291 ] . Certainly we have a slightly larger central value than achieved in 18j, where the same 
data are combined with the LCSR result f D f+^(fy = 1371$ MeV. 



To proceed, substitution of (56) in (38) gets 

f^(0) = 0.62 + 0.03, 



(57) 



where all the theoretical and experimental errors are in quadrature covered in the total 
uncertainty. There are abundant researches on "^(0), from which we pick just out several 
typical LCSR and lattice predictions and arrange them, along with the present estimate, into 
the tabulation in Tab. 2. At first sight, there exists a good accordance among all the LCSR 
results listed. Yet this should not be taken too seriously, for more or less parameters, on 
which the sum rules have relatively sensitive dependence, are chosen to have different inputs 
in these calculations. For instance, the obviously different parameter values are employed 
for both fu and 02 in jlTJ and 18] . The LQCD evaluations turn out to have a different 
extent of deviation from each other in the central values, but without any conflict within 
the errors. A comprehensive survey shows that f+~^ n (0) prefers taking a value larger than 
0.6. 

Lastly, we would like to present, as a by-product, our assessment for the decay constant 
fjj. From (50) and (56) follows that 



f d = 190+if MeV 



(58) 



with the same error disposal as in the f+~*"*(0) case. It falls into a somewhat wide interval 
formed by the existing findings of fo, which can be illuminated by the following examples. 
The CLEO measures f D = 205.8 ± 8.9 MeV fl on the assumption \V cd \ = \V US \ = 0.2255 ± 
0.0019. LQCD simulation predicts the three-flavor results f D = 218.9 ± 11.3 MeV 34] and 
f D = 213 + 4 MeV [35], and two-flavor one f D = 197 + 4 MeV Q. Compared with aUthese 



determinations, QCD sum rules provide, besides the two-loop result 203 ± 20 MeV (36j, the 



three-loop ones f D = 195 + 20 MeV [37| and f D = 177 + 21 MeV [38J. Hence one should step 
up efforts to improve calculations and promote understanding of that quantity. The present 
estimate, however, could be accommodated by fo ~ 200 MeV, a result gradually becoming 
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accepted on the basis of a multitude of phenomen olog ical investigations, and in particular 
accords well with those from the two-flavor LQCD 33J as well as three-loop QCD sum rules 



37j |. Meanwhile, these consistencies further expand support for the validity of our findings 



in (56) and (57). 

In the above discussion, a cross check has been made automatically between our LCSR 
predictions for /+~*' 7r (0) and /+~*' 7r (0). In contrast, it is out of the question for a traditional 
LCSR calculation, since twist-2 and -3 contributions, as emphasized, dominate respectively 
in the two sum rules, which consequently show exceedingly different sensitivities to both of 
them. In addition, from the observation that the twist-2 part predominates entirely over the 
twist-4 one in the present LCSR framework, we can benefit a lot in attempting to acquire 
a constraint on a 2 and 04 from the data on B and D decays. No doubt, this would enhance 
significantly our confidence in the LCSR applications to heavy-to-light transitions. 



IV. SUMMARY 

We have addressed in some detail a QCD assessment for B, D — > 7r transitions at the 
zero momentum transfer, in an improved LCSR approach, and presented our determinations 
of the form factors 7+^(0) and /^ )_i " r (0) as well as the CKM matrix elements \V u b\ and 
\Vcd\- We have also yielded a numerical estimate of the decay constant fjj. 

To 0(a s ) accuracy for twist-2 contributions and with the MS masses for the heavy 
quarks, the LCSR calculation on /bZ+^CO) and /£>/+ _ *' 7r (0) is carried out and the resulting 
sum rules bear the two remarkable characteristics: (1) They receive no contribution from 
not only the twist-3 but also the unknown twist-5 components, which are regarded usually 
as a serious source of uncertainty in the conventional LCSR applications, among others, to 
D decays, and therefore are available to twist-5 accuracy. (2) The twist-2 parts play a fully 
dominant role over the twist-4 ones so that the twist hierarchy required for convergence of 
the light cone expansions is preserved well and the higher-twist effects are kept under good 
control. The numerical analysis is performed with the updated inputs and experimental 
data; the validity and the self-consistency of the sum rule results are checked up and verified 
by a numerical comparison with some of typical theoretical predictions. Our findings are 
such as below: 

/B-ht(o) = 0.281°;°*, \V ub \ = (3.4+°^ ± 0.1 ± 0.1) x 10" 3 , 
/^(O) = 0.62 ± 0.03, \Vcd\ = 0.244 ± 0.005 ± 0.003 ± 0.008, 
f D = 190±i? MeV. 

The present results can be improved once the related inputs or experimental data become 
updated. Albeit unlikely to give help in understanding the existing discrepancy between in- 
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elusive and exclusive |T4&| determinations, an improvement on the \V u b\ determination is 
expected especially. However, it demands evidently a more decided knowledge of apart 
from a significant advance in experiment and in theoretical or phenomenological research 
on the pionic twist-2 DA. A continued and intensive study of fs helps also in the confir- 
mation whether or not non-SM physics shows an explicitly observable effect in r-leptonic 
and corresponding semileptonic B decays, which are expected to be detectable to a high 
precision in the running LHC or foreseeable super B factor. On the other hand, although so 
far our discussion on the form factors has been restricted to the largest recoil point q 2 = 0, 
q 2 dependence of them is understandable within the kinematical regions allowed by their 
individual light-cone expansion calculations. Then it is possible to extrapolate the results 
to the large q 2 regions in various ways available so as to have an all-around understanding 
of their behaviors. Too it is interesting to generalize the present discussion to the decays 
into K meson. We put off these studies to a future issue. 
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